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Recovery after cardiac surgery is a complex process that has to compensate for both individual variabiUty 
and extensive tissue damage in the context of systemic inflammation. Protein glycosylation is essential in 
many steps of the inflammatory cascade, but due to technological limitations the role of individual variation 
in glycosylation in systemic inflammation has not been addressed until now. We analysed composition of 
the total plasma and IgG N-gly comes in 107 patients undergoing cardiac surgery. In nearly all individuals 
plasma N-glycome underwent the same pattern of changes in the first 72 h, revealing a general mechanism 
of glycosylation changes. To the contrary, changes in the IgG glycome were very individualized. 
Bi-dustering analysis revealed the existence of four distinct patterns of changes. One of them, characterized 
by a rapid increase in galactosylated glycoforms, was associated with nearly double mortality risk measured 
by EuroSCORE II. Our results indicate that individual variation in IgG glycosylation changes during acute 
systemic inflammation associates with increased mortality risk and indicates new avenues for the 
development of personaUzed diagnostic and therapeutic approach. 

N- -glycosylation of human proteins is an essential posttranslational modification' generated by a complex 
biosynthetic pathway comprising hundreds of glycosyltransferases, glycosidases, transcriptional factors, 
ion channels and other proteins^. This process results in the creation of branched oligosaccharide chains, 
called glycans, which become integral part of proteins and significantly contribute to their structure and 
function"". 

Structural details of the attached glycans are of great physiological significance and many pathological con- 
ditions are associated with various types of glycan changes''"\ Glycosylation appears to be particularly important 
in the immune system' and immunoglobulin G (IgG) is one of the best-studied glycoproteins. Each Ch2 domain 
of the Fc region of IgG heavy chains carries a covalently attached bi-antennary N-glycan at the highly conserved 
asparagine 297 residue^. These glycans are essential for all pro-inflammatory activities of antibodies by main- 
taining the heavy chains in an open conformation required for FcyR interactions and even a small change in their 
composition has dramatic consequences for effector function of IgG" '. The addition of fucose to the core of these 
glycans radically reduces IgG binding to FcyRIIIa'"'" preventing initiation of antibody dependent cellular cyto- 
toxicity (ADCC)'^ ''' and destruction of target cells. On the other hand, the addition of sialic acid converts IgG 
from pro-inflammatory to anti- inflammatory agent" '^. Instead of binding to FcyRs, sialylated Fes initiate an 
anti- inflammatory cascade involving the lectin receptor DC-SIGN"" '^. This leads to upregulated surface express- 
ion of the inhibitory FcyRIIb on inflammatory cells, thus attenuating autoantibody-initiated inflammation. 
Sialylation of IgG was also found to be essential for the function of intravenous immunoglobulin (IVIG)'", a 
therapeutic preparation of highly purified polyclonal IgG antibodies widely used for the treatment of a number of 
autoimmune diseases. 

The importance of individual variation in glycosylation in acute inflammation has not been extensively 
studied, primarily due to the absence of reliable quantitative methods for high-throughput glycomics. We have 
recently developed methods for quantitative analysis of plasma" and IgG^" glycomes and here we have applied 
them to follow the dynamics of glycosylation changes during early course of systemic inflammation caused by 
cardiac surgery. Acute systemic inflammation is a part of many pathological events and patient's inflammatory 
response often determines outcome of a disease. Controlling the inflammatory cascade still represents a great 
challenge, due to its complexity and individual physiological differences. Cardiac surgery is a procedure that 
provokes a vigorous inflammatory response, which has important clinical implications^'"^'. This inflammatory 
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response is caused by blood contact with foreign surfaces and the 
activation of complement. Besides the activation of the complement 
system, increasing production of cytokines, oxygen radicals, release 
of endothelin and the expression of adhesion molecules on leuko- 
cytes and the endothelium represent main molecular mechanisms of 
such inflammation. Cardiovascular surgery with cardiopulmonary 
bypass (CPB) has improved in past decades, but inflammatory 
activation in this setting is still unpredictable and is associated with 
several postoperative complications^''. It is generally accepted that 
CPB initiates a whole-body inflammatory reaction. Factors influ- 
encing incidence, severity, and clinical outcome of the inflammatory 
response, and in particular the reasons why certain patients develop 
life-threatening perioperative complications, are currently not well 
understood. 

EUROSCORE (European System for Cardiac Operative Risk 
Evaluation) is a risk model which allows the calculation of the risk 
of death after a heart operation. The model asks for 17 items of 
information about the patient, the state of the heart and the proposed 
operation^"', and uses logistic regression to calculate the risk of 
death^"". This model has been adopted worldwide, becoming the most 
widely used risk index for cardiac surgery""^, and its use is believed to 
have contributed substantially to the improvement in the results of 
heart surgery seen at the beginning of the millennium. 

Due to enhanced production of cytokines in cardiac patients com- 
pared to other surgical procedures and modulatory effects of IVIG on 
the cytokine network, the supplemental application of immunoglo- 
bulins represents a potential therapeutic concept in cardiac post- 
operative high risk patients'". From the results of different studies 
it seems that immunoglobulins are unlikely to improve outcome in 
the entire group of patients, but only in specific subgroups. There is a 
possibility that dynamics and potential individuality of glycosylation 
changes analyzed in this work should be considered when the immu- 
nomodulation therapies are subjected. 

Results 

Glycosylation of plasma proteins. Total plasma N-glycome has 
been analysed in 107 individuals (Table 1) at three time-points: (1) 
prior to surgical procedure (day 0), (2) on the first postoperative day 
(day 1) and (3) on the third postoperative day (day 3). N-linked 
glycans were released from plasma samples using PNGase F and 
separated into 16 chromatographic peaks (Figure 1) as we have 
described previously". Since all patients have received heparin 
during the surgical procedure, the first six glycan peaks that 
showed aberrant chromatogram pattern and presumably contained 
degradation products from heparin have been excluded from the 
glycome for quantitative analysis. Since these chromatographic 
peaks contain glycans that nearly exclusively originate from IgG^*, 
and the IgG glycome was analysed separately, the exclusion of 
these peaks from the plasma glycome did not lead to any loss of 
information. 

Individual changes of plasma glycome composition in all patients 
are shown in Figure 2 (average values and their statistical significance 
are available as Supplementary Table 1). In nearly all analysed indi- 
viduals glycan groups GP9 (biantennary sialylated structures and 
oligomannose Man8 structure) and GP16 (tetra-antennary fucosy- 
lated structures) increased rapidly after the onset of systemic inflam- 
mation, while levels of aU other glycan groups decreased. Since we 
have previously demonstrated that individual plasma glycome has 
great temporal stability (average CV = 5.56%)^'^, the observed 
changes are apparently related to the disturbance of homeostasis 
induced by acute inflammation (reproducibility of the analytical 
procedure is CV = 7.4%)"'. 

Glycosylation of IgG. Contrary to the plasma glycome, temporal 
stability of the IgG glycome within an individual was not examined 
to date, thus prior to studying changes in IgG glycans after cardiac 



Table 1 Demographic and perioperative data 
N 


107 


Aqe (years) 




O/ 1 — 


Male (n/%) 




o 1 ^/ oj 


BSA (m^) 




9-1-0 9 


LVEF (%) 




— 1 o 


Smokinq (n/%) 




14 (131 


Hypertension (n/%) 






Diabetes (n/%) 






Chronic pulmonary disease (n/%) 






Renal failure (n/%) 




0 (/J 


Preoperative medicotion (n/%} 






Beta blocker 




74 (69) 


Statin 




70 (65) 


ACE-I 




67 (63) 


Operative procedure (n/%] 






Isolated valve surgery 




28 (26) 


Isolated coronary bypass surgery 




61 (57) 


Combined procedures* 




1 R 1 7 


Perioperative data 




Myocardial ischemia (min) 




78 34 


Cardiopulmonary bypass (min) 




112 ± 51 


Postoperative stroke (n/%) 




2(2) 


Perioperative myocardial infarction 


(n/%) 


3(3) 


Renal replacement therapy (n/%) 




3(3) 


30-day mortality (n/%) 




4(4) 


Highest post-operative CRP median 


(min-max) 


148 (36-399) 


BSA, body surface area; LVEF, left ventricular ejection fraction; ACE-I, angiotensin converting 
enzyme inhibitor. 

*lncludes multiple valve surgery, valve + coronary surgery, valve + aortic surgery. 



surgery, we evaluated "normal" day to day changes by analysing IgG 
glycosylation in 9 healthy individuals (age 22-78 years) through four 
days (day 1, day 2 and day 4). Glycans released from IgG molecules 
were separated into 24 chromatographic peaks (GP1-GP24), 
individual glycan structures were assigned to each peak and 
quantified as described previously^". A representative chromato- 
gram with glycan structures assigned to each chromatographic 
peak is shown in Figure 3. The IgG glycome within all individuals 
was remarkably stable, with average coefficient of variation of 5.9% 
(Table 2), what is comparable to variation in repeated measurement 
of the same sample. 

Changes in the IgG glycome composition after cardiac surgery 
were investigated by analysing N-glycans released from isolated 
IgG (and were thus not affected by changes in the concentration of 
IgG in plasma). Contrary to the plasma glycome which underwent 
similar changes in nearly all individuals (Figure 4A), changes in the 
IgG glycome were very individualized (Figure 4B and Supple- 
mentary Figure 1). Despite being discordant between different indi- 
viduals, these changes were surprisingly extensive and reached over 
100% change for many glycans and individuals (see Supplementary 
figure 1). In the context of significant intra-individual stability of the 
IgG glycome composition (Table 2, CV = 5.9%), this dramatic 
change suggests an important role of IgG glycosylation in the studied 
pathophysiological condition. 

Aiming to identify individuals with similar patterns of changes we 
have applied a bi-clustering algorithm, which in the same time clus- 
ters changes in individual glycans and groups individual patients 
according to similarity in patterns of glycan changes. Four stable 
groups (A - D) were identified using this algorithm (Figure 5A). 
The comparison of these groups with the available clinical data 
showed that one of the groups has a nearly double average 
EuroSCORE II (Figure 5B), a composite index used to predict risk 
of mortality after cardiac surgery'"'", that results significantly higher 
than the other groups (P = 0.0068, 1-way ANOVA performed with 
log(EuroSCORE) values, that have been used due to the nongaussian 
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Figure 1 | Part of the chromatogram of 2-AB labeled JV-linked glycans released from plasma proteins and separated by HILIC-HPLC. The integration 
areas, together with structures present in each glycan group are given. Groups of glycans are numbered from GP7-GP16. 



distribution of EuroSCORE values). Interestingly, some clinical vari- 
ables can be significantly associated with group 2: the number of 
patients that take diuretic drugs is significantly higher (P = 0.018, 
Kruskal-Wallis test), and also a composite index consisting in the 
number of drugs taken (among 'beta blockers', 'ACE inhibitors', 
'diuretics', 'statins', 'amiodaron') is significantly higher in group 2 
(P = 0.034, Kruskal-Wallis test). Closer examination of changes in 
IgG glycans in individual patients from this group (Figure 5C) 
revealed that they are characterized by a rapid decrease in IgG gly- 
cans without galactose (GP6: P = 0.0004, performed with a coupled 
Student's T test with a Post-Hoc Bonferroni correction for 96 tests 
over 24 peaks and 4 groups) and an increase in IgG glycans with two 
galactoses (GP14, P = 0.0086) and glycans with two galactoses and 
one sialic acid (GP18, P = 0.0067). These changes imply the more 
rapid increase in the antiinflammatory glycoforms of IgG (galacto- 
sylated'^ and sialylated^"'^') in group 2, than in other groups of 
patients. Moreover, considering the differences of the values for each 
peak between day 1 and day 0, the second peak (GP2) resulted 
significantly higher in group 2 as compared to the others (P = 
0.023, 1-way ANOVA, see Supplementary figure 2).When we tried 
to correlate log(EuroSCORE) with the glycan peaks at different times 
(Day 0, 1, 3) and with the differences (Dayl-DayO, Day3-Day0), 
some peaks show highly significant correlations/anticorrelations 
with the score, even if the values are not very high (|rj ~ = 0.4 for 
the highest correlations; Supplementary table 4.). 

N-GIycans as predictors of inflammation severity. When pheno- 
typic characteristics that indicate the rate of the inflammatory 
response (highest CRP during post-operative hospitalization, 
number of leukocytes and number of thrombocytes) have been 
compared to the preoperative levels of main characteristics of IgG 



glycans previously associated with inflammatory diseases (level of 
sialylation, fucosylation and presence of bisecting N-acetylglucosamine), 
statistically significant correlations for fucosylation were observed 
(Table 3). Statistically significant associations were also observed 
for the extent of fucosylation change between day 0 and day 1. 
Changes in the level of fucosylation can explain over 15% (R^ = 
0.15) of the CRP measured, what corresponds to correlation of 
approximately r = 0.4 (Table 3). These results indicate that, for 
example, a patient with high preoperative fucosylation level or 
moderate fucosylation change would have milder inflammatory 
response than a patient with low preoperative fucosylation or 
more prominent change in the level of this glycan feature. Since 
previously shown connection between glycans and age^" suggests 
that preoperative glycosylation could also be dependent on 
patient's age, change in fucosylation between day 0 and day 1 has 
higher clinical potential as prognostic marker. Since opportune 
intervention is of a great importance the day 0 and day 1 samples 
are taken into account for the creation of the prognostic marker, 
whilst the day 3 was neglected. 

Discussion 

Cardiac surgery provokes a vigorous inflammatory response with 
important clinical implications^'-^^, but factors influencing incidence, 
severity, and clinical outcome of the inflammatory response, and in 
particular the reason why certain patients develop life-threatening 
perioperative complications, are currently not well understood. Our 
recent studies revealed significant inter- individual variation of both 
total plasma (CV = 24.03%) and IgG glycome (CV = 39.4%) com- 
position''^'^". Since these values are much higher than those observed 
for method reproducibility, the inter-individual variation in IgG 
glycosylation is a consequence of the physiological differences within 
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Day 

Figure 2 | Changes in AT-linked plasma glycans through the initial phase of acute inflammatory response induced by cardiac surgery. Day 0 represents 
level of glycan group before the surgical procedure, day 1 on the first day after the surgical procedure and day 3 on the third day after the surgery. (The 
same patient is not represented by the same colour in all graphs, but a same spectrum is used for better visualisation.) 
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Figure 3 | Chromatogram of 2-AB labeled JV-linked glycans released from IgG molecules and separated by HILIC-UPLC. The integration areas, 
together with structures present in each glycan group are given. Groups of glycans are numbered from GP1-GP24, as used in the paper. 
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Table 2 | Intra-individual stability of IgG glycans, given as coefficients of variation of all HILIC UPLC separated glycan 
through four days in nine individuals 


groups measured 


Glycan 
group 












Coefficient of variation [%) 












Person 1 


Person 2 


Person 3 


Person 4 


Person 5 Person 6 


Person 7 


Person 8 


Person 9 


average 


reproducibility of 
the analytical 
inethod* 


GPl 


1 2.6 


9. 1 


1 2 


1 6.7 


20 


0.09 


0 1 

7 . 1 


O A 

0.4 


1 4.9 


1 1 .5 


Zo./ 


GP2 


5.06 


8.6 


8.9 


4.56 


3.6 


2.77 


o .o 


8.1 


5.22 


5.62 


o oi 

7 .7/ 


GP4 


1.16 


0.9 


3.2 


2.06 


2.6 


3.39 




1 . 1 


3.07 


2.34 


1 oo 


GPS 


7.01 


2.9 


2 1 


5.52 


8 


6.34 


o.o 


7.8 


8.3 


8.1 8 


o.o/ 


GP6 


0.94 


1 . 1 


4.2 


0.49 


3.4 


1 .6 


O . 7 


1 .5 


2.91 


2.23 


1 1 A 
1 1 .0 


GP7 


1 0.6 


2.6 


7.5 


4.28 


4 


4.38 


£; 1 
J . 1 


1 .2 


1 0.6 


5.59 


AAA 


GPS 


0.94 


0.8 


2 


0.66 


0.7 


0.74 


1 0 

1 .7 


1 . 1 


1 .09 


1 .09 


O ^A 
7 .DA 


GP9 


0.45 


U.4 


A "7 
4./ 


1 .78 


2.1 


2.01 


9 A 
z .o 


2.6 


3.48 


2.23 


T HA 


GPl 0 


2.35 


1 .3 


1 .6 


0. 1 8 


1.2 


1.75 




2.2 


0.39 


1 .54 


O AT 


GPl 1 


6. 1 4 


1 .6 


A A 
A.A 


4.59 


2.1 


0.98 




3.2 


6.2 1 


3.63 


0 TA 


GPl 2 


1 .2 1 


2.9 


1 3 


5.06 


5 


4.16 


U. 7 


8.1 


2.1 7 


4.73 


A OR 
4.V0 


GPl 3 


6.5 1 


5.8 


1 5 


0.97 


10 


18.3 


O.O 


9.6 


1 8.8 


9.84 


A 

O.OO 


GPl 4 


1.14 


0.5 


1 .3 


1 .33 


1 


0.81 


0.7 


1 .3 


1 .67 


1 .08 


/ .3(5 


GPl 5 


0.85 


0.5 


1 .5 


1 .56 


1.9 


2.15 


0.8 


1 . 1 


3.78 


1 .59 


T OA 


GPl 6 


1 .24 


0.6 


1 .1 


1 no 
1 .93 


3.5 


3.63 


5.3 


5.9 


1 .83 


oil 
3. 1 1 


A 00 
4.ZZ 


GPl 7 


4.05 


1.8 


4.7 


6.19 


4 


3.16 


6.7 


3.7 


8.15 


4.72 


T T3 


GPl 8 


1.05 


0.7 


3.5 


0.68 


1.9 


1.4 


3.3 


2.9 


3.66 


2.1 1 


A no 

O.UZ 


GPl 9 


2.45 


2.4 


5.1 


3.86 


3.1 


3.07 


7.5 


3.3 


7.95 


4.31 


A A"^ 
4.40 


GP20 


7.73 


5 


4.9 


15.6 


9.7 


10.1 


30 


7.3 


6.5 


10.7 


5.08 


GP21 


11.1 


8.9 


29 


13 


14 


8.25 


24 


13 


12.7 


14.8 


8.54 


GP22 


15.7 


10 


21 


4.56 


40 


12.1 


37 


21 


13.4 


19.4 


1 5.9 


GP23 


3.72 


7.4 


10 


0.74 


5.5 


3.6 


26 


8 


8.2 


8.18 


14.2 


GP24 


5.97 


8.1 


8.7 


0 


3.6 


5.58 


22 


12 


5.31 


7.9 


9.08 


average 


4.78 


3.76 


8.25 


4.18 


6.56 4.37 


9.12 


5.81 


6.53 


5.93 


8.65 


*determined by ten repeatec 


analysis of the some sample (IgG 


isolation, glyca 


1 release and labelling, HPLC 


analysis) on different occasions. 









A 




Figure 4 | Changes in levels of plasma glycans (A) and IgG glycans (B) between day 0 and day 1. Changes are given as log2 day 1/day 0. The red 
colour indicates negative values (decrease in glycan level) and blue colour indicates positive values (increase in glycan level); intensity of the colour 
indicates the intensity of the change. 
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Day 

Figure 5 | Bi-clustering analysis of patients and changes in the IgG glycome composition. (A) Bi-clustering heat map; (B) EuroSCORE II in 
subpopulations generated by bi-clustering; (C) individual changes in IgG glycans in subpopulations generated by bi-clustering. 



population rather than limitations of the method. Since glycans have 
essential structural and functional roles on numerous proteins'"', 
individual variations in glycosylation are expected to have significant 
functional implications^'^' This study is the first large study of 
individual changes in glycosylation of human IgG and total plasma 
proteins upon initiation of major inflammation. We used pre- 
inflammatory samples from each individual as a baseline control, 
thus eliminating the problem of the large inter- individual variability 



in glycosylation that burdens cross-sectional studies. Therefore even 
though our study was not very large (107 patients), each post-oper- 
ative sample was matched to a perfect control (pre-operative sample 
from the same individual). 

We observed extensive changes in the glycosylation of plasma 
proteins after cardiac surgery. Previously we have shown that com- 
position of the plasma glycome is characterized by temporal stabi- 
lity^'^, but it can apparently change very quickly if the homeostasis is 
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Table 3 | Connection of glycosylation features related to fucosylation, sialylation and bisecting GlcNAc with measures of severity of 
inflammatory response (highest CRP, number of leukocytes, number of trombocytes) 



Glycan feature 


Phenotype 


R 




P 


p adjusted 


roun rorai/oun 


vr ODD 

max (^Kr 


n TOT 
— U.oVO 


U. 1 JD 


<U.UU 1 


A AAt 


Fn total 


» ODD 

max LKr 


— U.oO 1 


A 1 T 1 
U. 1 O 1 


U.UUl 


A AAA 

O.OOO 


r\j \ n total/O 1 n 


» — ODD 

max LKr 


— U.OOO 


A 1 1 T 
U. 1 1 O 




A A1 A 

O.0 1 o 


CDC /CC t-.!— 1 

rbo total/ho total 


thrombocytes 


— U.Z04 


A A"7A 
U.U/U 


O.OIO 


0.042 


CQ /c 


» ODD 

max LKr 


— U.z// 


A r\~7~7 

U.U// 


A A1 A 


A Ayio 
0.04.Z 


rCj7zn total/Ozn 


max LKr 


— U.zzo 


A ACO 


0.004 


U. 1 1 V 


c 
o 


leukocytes 


(J. 1 vZ 


A AT7 


U.UO 1 


r\ 1 o 1 
U. 1 0 1 


c 
o 


max L-Kr 


U. 1 Oo 


A AOT 


U. 1 Jo 


AT 
U.4 1 J 


roUn rorai/oun 


thrombocytes 


U.UcSO 


A AAT 


U.4 1 Z 


U.C500 


CQ /Q 

ro/o 


thrombocytes 


U.U<54 


A AAT 


U.4ZO 


U.boo 


ArOUn total/OUn 


max LKr 


A TOA 

U.oVU 


A 1 

U. 1 OZ 


A AAA 

o.ooo 


A AAO 

o.ooo 


AFGln total/Gin 


max CRP 


0.270 


0.073 


0.011 


0.081 


AFS/S 


max CRP 


0.264 


0.070 


0.014 


0.081 


AFG2ntotal/G2n 


max CRP 


0.21 1 


0.045 


0.050 


0.224 


AS 


leukocytes 


-0.165 


0.027 


0.109 


0.392 


AFGln total/Gin 


thrombocytes 


0.151 


0.023 


0.148 


0.443 


AFBStotal/FStotol 


thrombocytes 


0.140 


0.020 


0.177 


0.456 


AFG2ntotal/G2n 


leukocytes 


-0.107 


0.012 


0.299 


0.539 


AFGOn total/GOn 


thrombocytes 


0.1 12 


0.013 


0.283 


0.539 


AS 


max CRP 


-0.125 


0.016 


0.248 


0.539 



FGOn/GOn - the percentage of fucosylation of agalactosylated structures without bisecting GicNAc; 

Fn total - the percentage of all fucosylated structures in total neutral glycan fraction; 

FGl n total/Gl n - the percentage of fucosylation of monogaloctosyloted structures; 

FBS totoi/FS total - ratio of all fucosylated sialyloted structures with and without bisecting GicNAc; 

FS/S - level of fucosylation of sialyloted structures; 

FG2n total/G2n - the percentage of fucosylation of digoloctosylated structures; 
S - the percentage of sialyloted structures; 

FGOn totol/GO - the percentage of fucosylation of agalactosylated structures; 
Aindicotes the difference between day 0 and day 1 . 



challenged. Changes in the total plasma protein glycosylation within 
the first three days after the onset of systemic inflammation followed 
the same pattern in nearly all patients, indicating the existence of a 
general physiological response. Differences in protein glycosylation 
between patients and controls have been previously reported in dif- 
ferent inflammatory disorders^*". Our previous studies revealed 
changes in glycosylation during the early course of sepsis and acute 
pancreatitis^'"^' and the same changes were observed after the acute 
inflammatory response has been induced in the animal model"*". 
However, the design of all previous studies prevented them from 
identifying very early changes, while in this study we clearly demon- 
strated very extensive glycosylation changes in the first 24 hours 
after the onset of inflammation. 

The decrease in all sialylated glycans, with the exception of bian- 
tennary disialylated structures and tetraantenary fucosylated struc- 
tures (which increased), occurred within the first 24 hours of 
inflammation and continued to change in the same direction during 
the next 48 hours. The notable exceptions from this pattern of 
changes were triantennary tri sialylated structures with fucose on 
the outer arm (GP14, also known as the sialyl-Lewis X antigen), 
which decreased in the first 24 hours, but then significantly increased 
in the following 48 hours. Though complex and somewhat unexpec- 
ted, the same pattern was observed in nearly all studied individuals 
(Figure 2), indicating that this is the real biological phenomenon. 
Sialyl-Lewis X antigen has often been associated with inflam- 
mation""''^ and its structural analogue present on leukocytes is 
responsible for leukocytes' interaction with selectin ligands on 
endothelial walls and the initiation of inflammation'" The increase 
in this structure on plasma proteins after a rapid decline in the first 
24 hours seems to be a reflection of a molecular mechanism that 
controls inflammation by competitive inhibition of selectin adhe- 
sion. However, this part of our study is limited by the fact that plasma 
glycome composition is affected by both changes in glycosylation 
and changes in the concentrations of different plasma glycoproteins, 
thus we cannot distinguish between these two mechanisms. 



Contrary to the total plasma glycome, changes in the IgG glycome 
did not follow a common pattern. Instead, as it is clearly visible in 
Figure 4B, IgG glycome composition appears to change nearly at 
random after the onset of systemic inflammation. However, the 
extent of these changes in the normally very tightly controlled IgG 
glycome'"''"' is evidently a consequence of specific physiological 
mechanisms resulting in significantly altered composition of the 
IgG glycome. Bi-clustering analysis of all patients revealed the exist- 
ence of four stable groups of patients with distinct patterns of IgG 
glycosylation changes in the first 24 hours following the cardiac 
surgery (Figure 5). One of these groups, characterized by a rapid 
increase in IgG galactosylation, was associated with nearly double 
EuroSCORE II index. The technology we used to analyse glycans is 
currently still not routinely used in the clinics, but it is foreseeable 
that in the near future changes in IgG glycosylation in the first 
24 hours after surgery could also be used to identify individuals with 
higher mortality risk. Furthermore, these individualized changes in 
IgG glycosylation could also point to molecular mechanisms behind 
the increased mortality risk. Increased IgG galactosylation has an 
anti-inflammatory role'. The aberrant activation of a wrong glyco- 
sylation mechanism (suppression, instead of promotion of inflam- 
mation) could be a part of the molecular mechanism leading to 
increased mortality risk, and thus a potential target for the develop- 
ment of new therapeutic approaches. 

Molecular events behind the observed changes in IgG glycosyla- 
tion are not known, but they could include (i) rapid production of 
new, alternatively glycosylated molecules of IgG, (ii) specific removal 
of some IgG glycoforms from circulation, or (iii) specific re-glyco- 
sylation in the circulation. IgG half life in plasma is approximately 26 
days and since its concentration is plasma is on average 15 g/L, it 
seems highly unlikely that a 24 hour period is sufficient for so extens- 
ive changes in the glycome through the production of new IgGs. 
Previous studies of transferrin glycosylation"*'"' showed its rapid 
desialylation during inflammation, which can be explained by the 
activity of specific plasma sialidase or enhanced protein elimination. 
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However, the complex glycosylation changes observed on IgG can- 
not be explained by simple exoglycosidase action. Selective removal 
of some IgG glycoforms through the interaction of IgG with its 
receptors and the aggregation of immunoglobulins is a possible 
mechanism, but another intriguing possibility is the extracellular 
re-glycosylation of IgG. For decades it is known that numerous gly- 
cosyltransferases are specifically released in the circulation, but it was 
widely believed that concentration of activated sugar nucleotides in 
plasma is not sufficient for efficient glycosylation. Recently it has 
been convincingly demonstrated that thrombocytes can release suf- 
ficient quantity of activated nucleotides for efficient re-glycosylation 
at the site of blood clotting^^. Extracellular sialylation of IgG has also 
been suggested recently^" and the possibility of protein specific extra- 
cellular modification of glycosylation should be considered as an 
option during inflammatory response. 

The observed individuality in the IgG glycosylation changes dur- 
ing inflammatory response is also potentially important for the 
administration of immunoglobulin therapy and could be one of 
the molecular mechanisms that could explain why this therapeutic 
approach is efficient in only certain subpopulations of patients. Only 
a subset of IgG glycoforms are active component of this therapy'^'^^ 
and highly individualized and extensive glycosylation response is 
something that has not been evaluated in this context so far. 

Methods 

Study population. The institutional etliics committee of the University Hospital 
Center Zagreb approved the study protocol. AU consecutive cardiac surgical patients 
operated at the University Hospital Center Zagreb from October 2010 to February 
2011 were screened for the study. Patients undergoing elective cardiac surgical 
procedures employing cardiopulmonary bypass met the eligibility criteria. Patients 
with a left ventricular ejection fraction <30%, emergency procedures, and those 
unable to provide consent were excluded from the study. All patients provided written 
informed consent. 108 patients were enrolled in the study, but full set of samples were 
collected for a total of 107 patients who were then included in the final analysis. 
Relevant demographic, clinical and perioperative data are summarized in Table 1 . 
These data were taken from the patients' medical records and were not collected 
specifically for this study. From all patients blood was drawn on that day, before 
surgery (day 0), and on the first (day 1) and third day (day 3) following surgery. This 
period has been chosen due to the essence of the study to observe the early changes 
during postoperative acute response and to identify potential prognostic 
glycosylation based markers. The blood was drawn on blood tubes containing EDTA, 
in the morning, between 7 and 8 a.m. The blood was centrifuged and the plasma was 
separated and stored at the — 20"C until analysis. The study was approved by the 
University hospital Zagreb Ethical Committee and performed in conformance to the 
ethical guidelines of the Declaration of Helsinki. 

Analysis ofN-linked plasma glycans. Release of N-glycans from plasma proteins. The 
5 ml of blood was drawn on the EDTA as an anticoagulant, centrifuged and the 
plasma was separated. Plasma proteins from 10 (J.1 of plasma were reduced by adding 
dithiotreitol (DTT) (Merck, Darmstadt, Germany) to final concentration of 0.5 M, 
followed by 100 mM iodoacetamide (lAA) (Sigma- Aldrich, St. Louis, MO, USA) 
alkylation. Denaturated glycoproteins were incorporated in SDS-polyacrylamide gel 
blocks to minimize the possibility of renaturation and glycoprotein loss. Gels were 
alternately washed with 100% acetonitrile (J.T. Baker, PhOlipsburg, NJ, USA) and 
20 mM NaHCOg (Merck, Darmstadt, Germany) buffer causing continuous 
dehydration and hydration of the gels, respectively, to remove all free detergents, 
reducing and alkylating agents. The enzyme, N-glycosidase-F (PNGaseF, Prozyme, 
Hayward, CA, USA) was absorbed by gel blocks and used to release the N-glycans 
from the proteins. Glycans were removed from gels by alternate washing with water 
and 100% acetonitrUe, collected in a 96-well collection plate (Waters, Milford, MA, 
USA) and dried in a vacuum centrifuge (Thermo Scientific, Waltham, MA, USA)^^. 

2-AB (2-aminobenzoamide) labeling of released plasma glycans. Dried samples were 
exposed to 20 |J.L of 1% (v/v) formic acid, dried again and finally fluorescently labeled 
by adding 5 |iL of 2-AB labeling solution [55 mg anthranilamide (Sigma- Aldrich, St. 
Louis, MO, USA) + 66 mg sodium cyanoborohydride (Sigma- Aldrich, St. Louis, 
MO, USA) + 330 [xl glacial acetic acid (Merck, Darmstadt, Germany) + 770 [il 
DMSO (Sigma- Aldrich, St. Louis, MO, USA)]. Mixture was agitated for 10 min, 
incubated for 30 min at 65 C, agitated again for 10 min, and incubated for a further 
90 min. 

Cleaning and elution of labeled N-glycans. In this final step, glycans were immobilized 
on chromatography paper (GE Healthcare) washed with 100% acetonitrile to remove 
the excess of 2-AB and finally collected by elution with 1 .5 mL of water. Samples were 
completely dried and redissolved in 75 \iL of water for further analysis. 



Hydrophilic Interaction Chromatography (HILIC)-HPLC of plasma glycans. Released 
glycans were subjected to hydrophilic interaction high performance liquid chro- 
matography (HILIC) on a 250 mm i.d. 5 mm particle packed TSKgel Amide 80 
column (Tosoh Bioscience, Stuttgart, Germany) at 30'^C with 50 mM formic acid 
adjusted to pH 4.4 with ammonia solution as solvent A and acetonitrile as solvent B. 
Separation method used linear gradient of 65-53% acetonitrile (v/v) at flow rate of 
0.8 ml/min in a 60 min analytical run. HPLCs were equipped with a Waters tem- 
perature control module and a Waters 2475 fluorescence detector set with excitation 
and emission wavelengths of 330 and 420 nm, respectively. The system was calibrated 
using an external standard of hydrolyzed and 2-AB-labeled glucose oligomers from 
which the retention times for the individual glycans were converted to glucose units 
(GU). Glycans were analyzed on the basis of their elution positions and measured in 
glucose units then compared to reference values in NIBRT's "GlycoBase v3.0" 
database available at http://glycobase.nibrt.ie) for structure assignment. 

Analysis of IgG glycans. Isolation of IgG from human plasma. The IgG was isolated 
from the same plasma samples using protein G monolithic plates as described 
previously^". Briefly, 50 |j.l of plasma was diluted lOx with PBS, applied to the protein 
G plate and instantly washed. IgGs were eluted with 1 inl of 0.1 M formic acid and 
neutralized with 1 M ammonium bicarbonate. This method was previously shown to 
produce the highly purified IgG'^ (verified by SDS-PAGE). 

N-Glycan Release. Isolated IgG samples were dried in a vacuum centrifuge. After 
drying, proteins were denatured with addition of 20 \iL 2% SDS (w/v) (Invitrogen, 
Carlsbad, CA, USA) by incubation at 60' C for 10 min. Subsequently, 10 \iL 4% 
Igepal-CA630 (Sigma -Aldrich, St. Louis, MO, USA) and 0.5 mU of PNGase F in 
10 )J.L 5 X PBS were added to the samples. The samples were incubated overnight at 
37' C for N-glycan release. 

2-AB labeling. The released N-glycans were subsequently labeled with 2-AB. The 
labeling mixture was freshly prepared by dissolving 2-AB (Sigma- Aldrich, St. 
Louis, MO, USA) in DMSO (Sigma- Aldrich, St. Louis, MO, USA) and glacial 
acetic acid (Merck, Darmstadt, Germany) mixture (17:3 (v/v)) to a final con- 
centration of 48 mg/mL. A volume of 25 ]iL of labeling mixture was added to each 
N-glycan release sample in the 96-well plate. Also, 25 |J.L of freshly prepared 
reducing agent solution (2-picoline borane (Sigma-Aldrich, St. Louis, MO, USA) 
in DMSO - concentration of 106.96 mg/ml) was added and the plate was sealed 
using adhesive tape. Mixing was achieved by shaking for 10 min, followed by 
2 hour incubation at 65 C. Samples (in a volume of 100 |iL) were brought to 80% 
ACN (v/v) by adding 400 (iL of ACN. 

Free label and reducing agent were removed from the samples using HILIC-SPE. 
An amount of 200 |j.L of a 0.1 g/mL suspension of microcrystalline cellulose (Merck, 
Darmstadt, Germany) in water was applied to each well of a 0.45 ]im GHP filter plate 
(Pall Corporation, Ann Arbor, MI, USA). Solvent was removed by application of 
vacuum using a vacuum manifold (Millipore Corporation, Billerica, MA, USA). All 
weUs were prewashed using 5 X 200 ]iL water, followed by equilibration using 3 X 
200 |J.L acetonitrile/ water (80 : 20, v/v). The samples were loaded to the wells. The 
wells were subsequently washed 5 X using 200 |j.L acetonitrile/ water (80 : 20, v/v). 

Glycans were eluted 2 X with 100 |J.L (instead of 200 |,iL) of water and combined 
eluates were not dried, but either analyzed immediately by UPLC or stored at — 20^ C 
until usage. 

Hydrophilic Interaction Chromatography (HILIC)-UPLC. Fluorescently labeled N- 
glycans were separated by hydrophilic interaction chromatography on a Waters 
Acquity UPLC instrument (Milford, MA, USA) consisting of a quaternary solvent 
manager, sample manager and a FLR fluorescence detector set with excitation and 
emission wavelengths of 330 and 420 nm, respectively. The instrument was under the 
control of Empower 2 software, build 2145 (Waters, Milford, MA, USA). Labeled N- 
glycans were separated on a Waters BEH Glycan chromatography column, 100 X 
2.1 mm i.d., 1.7 |im BEH particles, with 100 mM ammonium formate, pH 4.4, as 
solvent A and acetonitrile as solvent B. Separation method used linear gradient of 75- 
62% acetonitrile (v/v) at flow rate of 0.4 ml/min in a 25 min analytical run. Samples 
were maintained at 5 C before injection, and the separation temperature was 60''C. 
The system was calibrated using an external standard of hydrolyzed and 2-AB labeled 
glucose oligomers from which the retention times for the individual glycans were 
converted to glucose units. Data processing was performed using an automatic pro- 
cessing method with a traditional integration algorithm after which each chro- 
matogram was manually corrected to maintain the same intervals of integration for all 
the samples. The chromatograms obtained were all separated in the same manner into 
24 peaks and the amount of glycans in each peak was expressed as % of total integrated 
area. 

Statistical analysis. Batch correction was performed for each glycan group sepa- 
rately using linear mixed model (using programme language R). In this model, 
dependent variable was log transformed fraction of glycan group in total glycome. 
Age, gender, and time of sampling were described as fixed effects, and individual 
differences and technical sources of variation were described as random effects. To 
get corrected values for fractions of glycan groups, estimated batch effect (tech- 
nical source of variation) was subtracted from each log transformed glycan 
measurement. 

Bi-clustering analysis. Clustering analysis was performed on the difference between 
IgG glycosylation peak intensity values at day 1 (after surgery) and day 0 (before 
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surgery) with the "clustergram" function (Matlab, Mathworks USA) with "average 
hnkage" option and "correlation" distance both for rows {corresponding to each 
glycosylation peak) and columns {all glycosylation peaks for each sample) defined as 
one minus Pearson's correlation coefficient (equal to zero for complete correlation 
and equal to two for complete anticorrelation between rows/columns). 
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